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ABSTRACT 

Extremely flat and inverted radio spectra as observed in galactic nuclei and BL 
Lac sources are still a challenge for fast particle acceleration models. Continuous 
acceleration by electric fields in reconnection regions can result in almost constant 
particle distributions and thereby in inverted synchrotron spectra independent 
from the details of the injected spectrum. These spectra are calculated from a 
spatially averaged diffusion transport equation in momentum space that includes 
systematic momentum gains and losses as well as finite particle escape lifetimes. 

Subject headings: Galaxies: nuclei - BL Lacaertae objects: general - Radio 
continuum: galaxies - Acceleration of particles 

1. Introduction 

Flat and even inverted radio spectra, 1^ ~ z/^" with a < 0, are a common feature of the 
non-thermal emission of galactic nuclei. For example Sgr A*, the very center of the Milky 
Way, as well as the galactic center arc region (Lesch and Reich 1992; Duschl and Lesch 
1994) exhibit an inverted spectrum with a > —1/3. The same holds for the nucleus of M81 
(Renter and Lesch 1996). Slee et al. (1994) found that compact radio continuum cores in 
about 70 % of radio emitting elliptical and SO galaxies have a fiat or inverted spectrum with 
a median spectral index of a = —0.3. Observations of the Seyfert galaxy NGC 1068 indicate 
an inverted spectrum with index a = —0.3, as well (Muxlow et al. 1996). For active galactic 
nuclei (AGN) radio surveys by Kiihr et al. (1981) and Rayner et al. (2000) have accumulated 
a huge wealth of data about fiat radio spectra. Of special interests are the fiat radio spectra 
of the BL Lac sources Mrk 421 (Tosti et al. 1998) and Mrk 501 (Edwards et al. 2000), both 
also strong TeV-emitters. However, AGNs and BL Lacs exhibit relatively steep power-law 
spectra at frequencies higher than the radio regime. The complete electromagnetic spectrum 
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of these objects is a result of different energy loss and gain mechanisms. Consequently, the 
total spectra show flat and steep power-law spectra up to the radiation cut-off (Miller et 
al. 1996; Terzian et al. 1999). Flat spectrum radio sources are known to be variable over 
the entire accessible part of the electromagnetic spectrum down to the shortest timescales 
probed so far (Wagner 1998). 

Flat spectrum radio spectra have been attributed to self absorption. For the self- 
absorption peak frequency to extend over some decades specific geometrical properties in 
the source are required - this is well known as "cosmic conspiracy" (Marscher 1977; Cotton 
et al. 1980). This model depends on the source size as 1^ oc z/r, where r is the radius of a 
conical geometry. It implies that the linear size of the source is inversely proportional to the 
emission frequency, which is not supported by interferometric data (Cotton et al. 1980). Of 
special interest is the galactic center region, in which it was clearly proven that its flat radio 
spectrum is not due to self absorption (Lesch et al. 1988; Yusef-Zadeh 1989). 

The evidence for flat, optically thin radio spectra in several active galactic nuclei has 
been presented by Hughes et al. (1989); Melrose (1996); Wang et al. (1997). All these 
authors consider different Fermi-like acceleration schemes (either multiple shocks or second 
order Fermi mechanism) to be responsible for the hardness of the electron energy spectra. 

In the present contribution it is suggested that optically thin synchrotron emission due 
to hard electron spectra produced in magnetic reconnection regions can explain the nature 
of the flat/inverted spectrum radio sources. 

2. Pcirticle acceleration in reconnection regions 

Any acceleration model of charged particles must ultimately be based on the energy gain 
of the particles in electric fields. However, the dynamics of the respective plasma processes 
involved can be quite different. The principal acceleration scenarios (reviews are given, e.g., 
by Schlickeiser (1986); Kirk (1994); Kuijpers (1996)) may be divided in the Fermi I and Fermi 
II mechanisms, shock-drift acceleration (which in fact can be considered as an example of 
Fermi I acceleration), diffusive shock acceleration, plasma wave acceleration, electrostatic 
double layers and the energization of particles in reconnection events. In contrast to the 
other 'classic' mechanisms it was only more recently that reconnection was discussed as a 
process that does not only convert magnetic field energy to plasma bulk motion and heating 
but also plays an important role in the context of fast particle acceleration (Schindler et al. 
1991; Vekstein and Priest 1995; Lesch and Birk 1998; Schopper et al. 1998; Litvinenko 1997, 
1999). 
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In the context of shock acceleration much elaborated work has been done concerning 
the shape and evolution of the energy distribution of the accelerated particles (e.g. Kuijpcrs 
(1996); Achterberg (2000) and ref. therein). One may note that in the diffusive shock 
scenarios the injection problem appears, i.e. one has to assume a high-energy electrons to 
start with. Within many applications only the modification of power law spectra by shocks 
which are assumed for the injected particle population in the first place are investigated. 
What is more, very flat spectra are very hard to explain in the context of simple shock 
acceleration (Melrose 1996). 

Nothing comparable to the detailed studies on the energy spectra of shock accelerated 
particles has been done for the reconnection scenario. In this contribution we analyze the 
energy spectra that can be expected to be caused by high-energy particles accelerated in dc 
electric flelds in reconnection regions starting from the space-averaged transport equation for 
the one-particle distribution function that represents the high-energy particle component. 

In three-dimensional configurations with B ^ everywhere a finite electric field com- 
ponent parallel to the magnetic field Eg is necessary for the onset of magnetic reconnection 
(Hesse and Schindler 1988; Schindler et al. 1991). Ideal Ohm's law is to be violated 



where R is some yet unspecified non-idealness. In the associated generalized electric potential 
(Schindler et al. 1988, 1991) 



where the integral is evaluated along the magnetic field lines that penetrate the reconnection 
region, charged particles can be efficiently accelerated. It should be mentioned that in general 
U is not an electrostatic potential. Note that particles are not necessarily accelerated up 
to energies associated with U but expression (2) gives an upper limit of the energy level 
of the accelerated particle population. The electric potentials are often associated with 
non-ideal flows like steady rotational motion or (non-)stationary shear flows. The nature 
of the non-idealness R and the reconnection process determines the strength, structure 
and time dependence of the electric potential and thereby the particle energization. One 
can think of different kinds of non-idealness in Ohm's law that allow for the reconnection 
process to start. Since astrophysical plasmas are on large spatial scales highly coUisionless 
the Coulomb resistivity is negligible. Nevertheless, Ohm's law can locally be violated by 
Rs — Tjjs, where jg is the magnetic fleld aligned component of the current density, provided 
that some anomalous resistivity rj is caused by microturbulent stochastic electromagnetic 
flelds (e.g. Huba (1985) and ref. therein). The microturbulence can be caused by currents. 



E + -vxB = Ry^O 



c 



(1) 




(2) 
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When the current density exceeds a critical value, unstable kinetic waves can be excited 
that in their non-linear evolution result in stochastic microturbulent electromagnetic fields. 
The same concept Rg ~ jg for j > jcru holds for stochastic ion acoustic double layers as 
a source for anomalous dissipation (Haerendel 1993). In relativistic plasmas, in particular, 
finite particle inertia can be responsible for the non-idealness Rg ~ djs/dt + V ■ (vj + jv) as 
was recently discussed in the context of extragalactic jets (Lesch and Birk 1998). 

In the present study we are rather interested in the resulting power-law indices of the 
energy spectra of the high-energy particles than in their maximum energy. The spectra result 
from the combined effects of direct acceleration of charged particles in the dc electric fields 
in reconnection regions and radiative loss processes. For the acceleration the yet unspecified 
quantity Rg plays the crucial role. In the simplest standard scenario of two-dimensional 
stationary reconnection a constant electric field perpendicular to the plane where a projection 
of one magnetic field component changes sign occurs, e.g., (Parker 1963; Biskamp 1993). 
This electric field is a result of the non-idealness Rg that is assumed to be constant in the 
reconnection region. Let us now consider a configuration with a dominant constant magnetic 
field component parallel to the direction of the associated electric current. The electric field 
parallel to this magnetic field component is approximately Eg and therefore the fast particle 
acceleration is constant in time. A finite pitch angle of the accelerated particles allows for 
the radiative losses. To be more specific let us consider a generic shear configuration. If a 
plasma fiow, say = Vg{x,z)ey is applied to an homogeneous magnetic field B = Bzez a 
shear component forms By — J Bzdvg/dzdt that results, in particular, in an electric current 
fiowing parallel to B^. Given some localized dissipation (caused, e.g., by a supercritical 
current density) magnetic reconnection can take place. In a stationary state a constant 
(Biskamp 1993) inflow Vj = ViG^ towards the reconnection region occurs and the magnetic 
shear is converted to heat and kinetic energy as fast as it is built up by the external shear 
flow Vg. The convergent velocity Vi can be expressed as a function of the local Alfvcn velocity 
and the (constant) reconnection rate (Biskamp 1993). The quasi-stationary magnitude of the 
electric field component that causes the acceleration of fast particles along Bz can then be 
estimated without specifying for the non-idealness Rg as it is possible in a more elaborated 
general three-dimensional reconnection model (Schindler et al. 1988, 1991). In our example, 
we find 



where Arec and T denote the extent of the reconnection region along B^ and the typical time 
scale of shear formation. 

In what follows, we describe the particle acceleration by a an average constant electric 
field. In particular, the acceleration force does not depend on the momentum of the charged 
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particles. 



3. The momentum balance equation for the accelerated particle population 

The population of high-energy particles that do not interact with each other can be 
described by the one-body distribution function /(r, p, t) that obeys the Liouville equation 

Though this equation contains the full dynamical information it is of limited use since this 
non-linear differential equation in seven variables is usually unsolvable for physical situations 
of interest. Fortunately, we need not to know the full information of the particle dynamics 
in a rigorous manner. First, will assume pitch angle isotropy. This idealization is motivated 
by the fact that in reconnection regions Alfvenic and magnetosonic wave fields with rela- 
tively high-energy densities are to be expected which cause an efficient particle scattering in 
momentum space (Schlickeiser 1986). Additionally, self-excited Alfven turbulence maintain 
quasi-isotropy, if the systematic velocity of the high-energy particles exceeds the Alfven speed 
of the ambient plasma (Kuijpers 1996). Second, when asking about the energy spectrum of 
the accelerated particles in a specified astrophysical object the detailed spatial dependence 
of the distribution function is often of no interest. Rather it is usually sufficient to con- 
centrate on a distribution function space- averaged over the volume V. Thus, we deal with 
the isotropic distribution f'{p,t) = J drf(r,p,t)/V. The total particle number is calculated 
as N{t) — An J f'p^dp. The loss in information is not too high a price to pay in order to 
obtain a transport equation in momentum space one can really deal with. The continuous 
momentum gains and losses Ap the particles experience can be described by a momentum 
operator 

1 d 

fdp ^ Pioss)p'^] , (5) 

if one concentrates, as in the present study, on systematic temporal momentum gains Pgain 
and losses pioss- The finite characteristic time r a particle stays in the acceleration region 
gives rise to a further sink in the balance equation for f'ip,t). In the reconnection scenario 
particles can escape Eg ^ 0-regions by following reconnected field-lines outward the non- 
ideal regions. Some explicit source term S(p), on the other hand, represents the injected 
particle population. Taking these effects into account the temporal evolution of the particle 
distribution is governed by the equation (Jones 1970; Schlickeiser 1984, 1986) 

^ + £p/V*) + ^ = SW (6) 
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which holds, if the distribution function /(r, p, t) is separable in the space and the momentum 
variables, and if the momentum operator Cp is space independent. 

It should be noted that equations(5) and (6) hold only for momentum changes Ap <^ p 
during a time interval At ^ r. If one wants to include the effects of stochastic diffusion 
in momentum space, which may become important, e.g., for a microturbulent resistivity, a 
Fokker-Planck diffusion term could be incorporated (Schlickeiser 1984, 1986) in equation(5). 

The momentum change terms in equation(5) in our case can be specified as follows. The 
dominant loss processes that are mainly responsible for the non-thermal radiative phenom- 
ena are synchrotron and inverse Compton radiation. We will not consider Bremsstrahlung 
{pioss ~ p) here, since this effect usually is dominated by the other ones mentioned. Both 
processes scale as (Longair 1981) pioss = —kp^ with k = ^aTC^U where is the Thom- 
son cross-section and U is the magnetic field or photon energy density, respectively. The 
momentum gain of a single charged particle in the acceleration region is given by 

Pgain^ ^^{irnv) ^ qEs^qRs (7) 

where 7 is the Lorentz factor and q is the electrical charge. In this contribution we exclusively 
consider the case qRg — C =const. (see discussion in section 2) and thus, solve the transport 
equation 

|^+(^-«p)/+(C-;.p^)| + ^ = sw (8) 

where the primes are now omitted. Equations of this kind can, in general, be solved by 
means of the inverse Laplace transformation (Melrose 1980). 

4. Power law peirticle distributions 

First, wc solve the steady-state version of equation(8) for different source populations. 
The resulting distribution function reads 



f{p) = hp-\-kp' 

where /o depends on the injected source population by the integral 



/p 



l+pVk/C 



1-PvWC 



(9) 



p'^/kK 



S ip') dp'. 



(10) 
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If the acceleration is negligible as compared to radiative losses, the stationary energy spec- 
trum equation(9) shows the pure aging behavior (see Longair (1981)) caused by synchrotron 
or inverse Compton losses f{p) ~ where f{p)dp is the particle number in the momentum 
interval dp. If the losses are negligible, the spectrum is given by /(p) ~ p~^e~P/'^'^. 

We note that the distribution function given by equation (9) may diverge at the momen- 
tum p — Pc— -x/CA where gains and losses balance. When it is rewritten as 



\Pc 



1 + — 

Pc 



Pc 



(11) 



= 0. A pile- 
/2 if one defines 

-p/tloss = -kp"^, 



it becomes clear that the function is well behaved for 2Ty/k( < 1, with f{p, 
up of particles occurs if 2T\^k( > 1. This is equivalent to r > \/taccti 
the acceleration and loss timescale as Pgain = p/tacc = C Pioss = —p/i-ioss 
respectively. This means that the particles are collected near p = Pc when the escape time 
(which depends on the geometry of the acceleration region) is longer than half the geometric 
mean of the timescales for acceleration and radiative losses. Otherwise the particles escape 
before they can accumulate. 

It needs mentioning that the pile-up at the critical momentum p^ = \/C/k would be 
smoothed out, if a Fokker-Planck diffusion term was incorporated in equation(8). This was 
done by Schlickeiser (1984) for the case of shock acceleration. 

We now turn to determine /o in equation(9) for different injection spectra. For a mo- 
noenergetic source population S{p) = so6{p — po) we find 



/o = soplQ{p-po) 



'l+povwc^ 

i-poVWc 



(12) 



where denotes the Heaviside function. In the next section we will calculate spectra for 
this source population. For a relativistic Maxwellian source population S{p) = Soexp(— ap), 
where a is a constant, one obtains 



/o = So 



/' 



l-p'\/k/C 



dp' 



(13) 



where c is some integration constant. In the limit of negligible particle escape, i.e. Ty/k( 
oo, or for negligible acceleration as compared to radiative losses k/( ^ oo the particle 
distribution reads 

/(p) = -e-'^f (2 + 2ap + aV)|C - A;p=^r V- (14) 
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For a power law source distribution, S{p) — sop ^ fo, can be determined as 



fo = so r-'e 



i2-z -ap' 



i-p'aA7C 



dp'. 



(15) 



If one considers again the limits of negligible particle escape or negligible acceleration one 
finds 



■i-z 



fiP)^so^\c-kpr'p-'- 

6 — Z 



(16) 



Non-stationary solutions of equation(8) with an exponential time evolution /(p, t) ~ 
e"*/"^ can be found provided the distribution function satisfies Cpf — S{p). For a monoener- 
getic source population one obtains for p > po 



f{p,t) = cp \C-kp \ e 



(17) 



where the constant c contains Sq- A relativistic Maxwellian source population and a power- 
law one result in 

f{p,t)=(c+'-^e-'^Ap-'\C-kp'\-'e-^ (18) 

and 

f{p,t) ^ {c + Sop-') p-^\C - kp^\-'e-^ : (19) 

respectively. 



5. Synchrotron emission spectrum 

The spectral power P,y radiated by a single particle with momentum p is in a good 
approximation given by (b = y/Se^B sin 9/ mc^ and Uc = 3eB smOp'^/ATT^mc)^, with 9 being 
the angle between the magnetic field and the hne of sight) 

p.{p)^b(^^y\-''/''% (20) 

see Melrose (1980b). It has its maximum near u — Uf. and shows a characteristic u^^^ 
spectrum for u <^ i/c and an exponential cutoff for u ^ Uc- The specific intensity from an 
isotropic momentum distribution of particle density N{p) — 47rp^/(p) can then be calculated 

by 

4 = y" dpN{p)P,{p). (21) 
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In figure (1) we plot the distribution function N{p) — 47rp^/(p) with /(p) from equation(ll) 
for a monoenergetic source population (equation(12)) with different values of a = l/2T\fkC,. 
The pile-up results for a = 0.5. The higher cut-off results for a = 2 and the third curve 
represents the case a = 6. Figure (2) shows the corresponding synchrotron radiation spectra 
along with the spectrum of a single particle (upper curve) according to equation(20) for 
comparison. Going downward, the other curves again represent a = 0.5, a = 2, and a = 6, 
respectively. It is interesting that the low-frequency part becomes flattened even though 
the particle distribution remains constant at low energies and changes only near the cutoff 
momentum. The spectral index a with Ii, oc is given by a = —0.33 for the upper curve 
and a —0.1 for the lowest curve. This shows that synchrotron emission with inverted 
spectra is to be expected from reconnection zones. We note again that in the context of 
AGNs and BL Lacs the complete emission spectra show power-laws for higher than radio 
frequencies. 



6. Discussion 

Particle acceleration by field-aligned electric fields in reconnection regions can be de- 
scribed by diffusion in momentum space as long as the differential momentum changes are 
small compared to the respective instant particle momentum. This assumption seems not to 
be that restrictive, since reconnection events are expected to excite intense Alfvenic modes 
and MHD turbulence that causes efficiently pitch angle diffusion. Thus, in many applications 
the diffusion treatment seems to be appropriate. The present analysis shows that mainly 
constant distributions for accelerated particles are to be expected no matter whether mo- 
noenergetic, Maxwellian or power law populations are injected into the reconnection regions. 
The associate synchrotron spectra Ii, oc i/"" show negative spectral indices a < 0. Such ris- 
ing spectra are ubiquitously observed in galactic nuclei (Rayner et al. 2000) as in the core 
regions of blazars (Edwards et al. 2000). 

We stress that the actual acceleration problem is highly involved. Our findings come 
from a quite simple nonetheless physically reasonable scenario. Highly anisotropic particle 
distributions cannot be described within our approach. Neither local variations of the spectra 
can be determined. Whereas some improvement may be possible on the analytical side 
we feel that a numerical test particle approach and self-consistent electromagnetic particle 
simulations are indispensable for a deeper understanding of the acceleration of high-energy 
particles in reconnection regions. 
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Fig. 1. — The particle spectrum in a double-logarithmic representation for a — l/2Ty/k( (see 
equation(ll) chosen as a = 0.5 (pile-up), a — 2 (higher cut-off), and a — 6 (lower cut-off), 
respectively. 

Fig. 2. — The synchrotron emission spectrum in a double- logarithmic representation. The 
uppermost curve represents the single-particle spectrum and the following ones represent the 
cases a — 0.5, a — 2, and a — 6, respectively. The logarithm of the intensity is normalized 
to its value at log{u/i/c) — —6. 
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